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The analysis of nonperturbative effects in high energy asymptotics of the electomagnetic quark 
form factor is presented. It is shown that the nonperturbative effects determine the initial value for 
the perturbative evolution of the quark form factor and find their general structure with respect to 
the high energy asymptotics. Within the Wilson integral formalism which is natural for investigation 
of the soft, IR sensitive, part of the factorized form factor, the structure of the instanton induced 
effects in the evolution equation is discussed. It is demonstrated that the instanton contributions 
result in the finite renormalization of the subleading perturbative result and numerically are char- 
acterized by small factor reflecting the diluteness of the QCD vacuum within the instanton liquid 
model. The relevance of the IR renormalon induced effects in high energy asymptotic behaviour is 
discussed. The consequences of the various analytization procedures of the strong coupling constant 
in the IR domain are considered. 
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I. INTRODUCTION 



£> | The electromagnetic (color singlet) quark form factor is one of the simplest and convenient objects for the investiga- 

■ tion of the double logarithmic behaviour of the QCD amplitudes in the high energy regime. From the methodological 
point of view, the consistent study of such an asymptotic requires a perturbative resummation procedure beyond 
the standard renormalization group techniques. Besides this, the resummation methods developed for this particular 
case can be applied to study of many other processes which possess the logarithmic enhancements near the kinematic 
boundaries. On the other hand, in addition to the evident theoretical significance, the computation of the quark form 

■ factors has important phenomenological applications. The quark form factor enters into the cross sections of a number 
of the high energy hadronic processes [1]. For example, the total cross section of the Drell-Yan process (normalized to 
DIS) is determined by the ratio of the time-like and space-like form factors [2, 3]. The similar resummation approach 
is used also in study of the near- forward quark-quark scattering, and evaluation of the soft Pomeron properties [4]. In 
the latter case, the nonleading logarithmic terms are quite important. The investigation of the electomagnetic quark 
form factors (Dirac as well as Pauli) in moderate and low energy domains can shed a light on the problem of the 
scaling violation in DIS and the structure of constituent quarks [5]. 

The first example of the large logarithm resummation was given by Sudakov for the case of o/f-shell fermion in 
external Abelian gauge field in the leading logarithmic approximation (LLA), where the terms of order of (a™ In 2 ™ q 2 ) 
are taken into account while the contributions from 0(a™ In 2 " -1 q 2 ) are neglected. The exponentiation of the leading 
double logarithmic result was found [6]. This exponential decreasing of the form factor at large-*? 2 means that the 
elastic scattering of a quark by a virtual photon is suppressed at asymptotically large momentum transfer. The 
exponentiation for the cm-shell form factor in the Abelian case was obtained in the LLA in [7]. As expected, the non- 
Abelian gauge theories appeared to be more complicated: first, the leading LLA terms in the QCD perturbative series 
were found to be consistent with exponentiation in [8] (the inelastic on-shell form factor with emission of one and two 
gluons was calculated in the same context in [9]; the role of the quark Sudakov form factor in the description of e + e~ 
one-photon annihilation in quarks and gluons was considered in LLA in [10]), and the all-order LLA non- Abelian 
exponentiation has been proved in [11]. In the LLA, the exponentiated form factor was shown to be the rapidly 
decreasing function at high momentum transfer, but the question if the non-leading logarithmic terms could upset 
the LLA behaviour required a further work. The all-(logarithmic)-order resummation was performed in the Abelian 
case and the exponentiation was demonstrated in [12]. In the paper [13], the non- Abelian all-order exponentiation for 
the so-called hard part of the on-shell form factor has been shown first within the powerful factorization approach. 
Note, that in this work the case where a time-like photon with large invariant mass decays into a quark-antiquark 
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pair was considered, however it can be easily shown that the results remain true for our case of a quark scattering in 
an external EM field as well. 

In the work [13], the detailed study of hard part of the form factor (which is responsible for the UV properties) was 
performed, while the status of the soft part, containing all the IR and collinear singularities and, as a consequence, 
all possible nonperturbative effects, remained unclear. The important results on the IR properties of the QCD vertex 
functions was obtained in [15] within the Wilson integral approach. In these works, the soft part of the form factor 
had been presented as the vacuum averaged ordered exponent of the path integral of a gauge field over the contour 
of a special form — an angle with sides of infinite length. The use of the gauge and renormalization group invariance 
allowed to derive the perturbative evolution equation describing the high energy behaviour of the form factor taking 
into account all (not power suppressed) parts of the factorized amplitude, both for the on- [14] and q/f-shell [16] cases. 
It was shown that the leading asymptotic is controlled by the cusp anomalous dimension which arises due to the 
multiplicative renormalization of the soft part, and can be calculated within the Wilson integral formalism up to the 
two-loop order [15]. It is worth noting that within the Wilson integral approach, the non-Abelian exponentiation can 
be proved independently [17], what is another important advantage of this framework. The efficiency of the Wilson 
integrals approach had been successfully demonstrated in a series of works [18-20]. In these papers, the consistent 
non-diagrammatic framework is developed what allows to calculate the fermionic Green's functions, Sudakov form 
factors, amplitudes and cross sections in QED and QCD completely in terms of world-line integrals, and thus avoid 
complicated diagrammatic factorization analysis. 

The results presented above allow one to conclude that the leading high energy behaviour of the quark form factor 
in non-Abelian gauge theory is completely determined by the perturbative evolution equation, and is given by the 
fast decreasing exponent: 

(1) 

This rapid fall off is not changed by any other logarithmic contributions [13, 14, 16, 19]. However, the non-leading 
logarithmic corrections are nevertheless important for evaluation of the numerical value of the form factor. Some of 
them are of a purely perturbative origin (higher loop corrections and sub-leading logarithmic terms), while the others 
can be attributed to the nonperturbative phenomena. The usual approach to treatment of the latter is developed 
within the IR renormalon picture (there are plenty of papers on this subject, for the most recent reviews see [21]). 
However they could only give the power-suppressed terms, which become, of course, important in low energy domain, 
but can be neglected at asymptotically large momenta. Here we should note, that this conclusion is to be changed for 
processes with two scales (such as quark-quark scattering, Drell-Yan process, etc.): then the corrections proportional 
to the powers of a smaller scale must also be involved in the game [22]. In the present work, we try to advocate the 
point of view that the true (not connected directly to renormalons) nonperturbative effects can be taken into account 
consistently in the evolution equation, and therefore they yield the non-vanishing subleading (perhaps, perimetrically 
suppressed, but still logarithmic) contributions <~ In q 2 to the high energy behaviour. Further, we analyze another 
possible source of contributions which can be considered as "nonperturbative" — the IR renormalon ambiguities. We 
demonstrate explicitly that they produce the corrections with different IR structure compared to that one generated by 
instantons. Moreover, as it can be shown these direct renormalon effects disappear in the dimensional rcgularization 
[3] and in the analytical perturbation theory [23] , what means that one could hardly expect a significant signature of 
the IR renormalon effects in this process. 

The idea that the nontrivial vacuum structure could be relevant in high energy hadronic processes was first ex- 
plicitly formulated for the soft Pomeron case in Abelian gauge theory by Low [24] , Nussinov [25] , and Landshoff and 
Nachtmann [26], and developed further using the eikonal approximation and the Wilson integral formalism in [27]. In 
the present work, we consider the nonperturbative effects originating in the nontrivial structure of the QCD vacuum 
treating the latter within the framework of the instanton liquid model (ILM) [28-31]. The approach based on the 
other principles is successfully developed within the stochastic vacuum model (SVM), where some important and 
interesting results have been obtained (see, e.g., [32]). However, since the correspondence between both pictures is 
not completely clear at the moment, we will restrict ourselves with the ILM only omitting the discussions of relations 
with results of SVM. 

The paper is organized as follows: In Section II we describe the consequences of the RG invariance properties 
of the factorized form factor, and derive the linear evolution equation considering the nonperturbative input as 
the initial value for perturbative evolution. In Section III, the nonperturbative effects are estimated in the weak- 
field approximation within the instanton model of the QCD vacuum. The large-g 2 behaviour of the form factor is 
analyzed taking into account the leading perturbative and instanton induced contributions. In Section IV, we study 
the consequences of the IR renormalon ambiguities of the perturbative series and discuss their relevance within the 
context of certain analytization procedures. 
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II. EVOLUTION EQUATION AND NONPERTURBATIVE EFFECTS 



The behaviour of the form factors in various energy domains is one of the most important questions in the theory 
of hadronic exclusive processes. The electromagnetic quark form factors are determined via the elastic scattering 
amplitude of a quark in an external color singlet gauge field: 

= F q [(pi -p 2 ) 2 } u( Pl )^v(p 2 ) - G q [(p! -p 2 ) 2 } u( Pl f^ iPl 2m P2) » v(p 2 ) , (2) 

where u(pi), v{p 2 ) are the spinors of outgoing and incoming quarks, and a^ v = [7^, 7„]/2. In the high energy regime, 
the Pauli form factor G q is power suppressed and will be neglected in the present consideration. However, it should 
be emphasized that in low and moderate energy domains it becomes important and interesting perturbative as well 
nonpcrturbative effects arise (see, e.g., recent works [33, 34]). 

The kinematics of the process can described in terms of the scattering angle x : 

cosh x = ^l = l + ^ , Q 2 = -(p 2 - Pl ) 2 >0 , p?=p§ = m 2 . (3) 

The classification of the diagrams with respect to the momenta carried by their internal lines allows to express the 
form factor F q in the amplitude (2) in the factorized form [12-14] (compare with the world-line expression for the 
three-point vertex in [19]) 

F q (q 2 ) = F H (q 2 /^ 2 , a s ) ■ F s (q 2 /m 2 , M 2 /A 2 , a.) ■ Fjtf/X 2 , a.) , (4) 

where the hard, soft, and collinear (jet) part are separated. Note, that in the present paper, all the dimensional 
variables are assumed to be expressed in units of the QCD scale A, so that q 2 = Q 2 /A 2 , etc. The arbitrary scale /s 2 
stands for the boundary value of the squared internal momenta which divides the different parts, and is assumed to 
be equal to the UV normalization point. 

The total form factor F q is the renormalization invariant quantity: 

^w)) = (^ 2 ^+^)^y ^(mW))=o , (5) 

what leads, in the large-g 2 regime, to the following relations 

! 



2 d 



' d In Fu ' 
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d In q 2 




d In q 2 



2 r cusp (a s ) . (6) 



For convenience, we work with the logarithmic derivatives in q 2 what allows us to avoid the problems with additional 
light-cone singularities in the soft part [14, 35]. The collinear part Fj being independent on q 2 does not contribute 
to these equations. 

Within the eikonal approximation, the resummation of all logarithmic terms coming from the soft gluon subprocesses 
allows us to express Fs in terms of the vacuum average of the gauge invariant path ordered Wilson integral [36] 

F s (q 2 /m 2 ,fj 2 /X 2 ,a s ) ^ W(C X ^ 2 /X 2 ,a s ) = ^Tr(0|7>exp ji ff jf dx^^j |0> . (7) 

In Eq. (7) the integration path corresponding to considering process goes along the closed contour C x : the angle 
(cusp) with infinite sides. The gauge field A^x) = T a A*(x) , Ti[T a T h ] = ±5 ab , belongs to the Lie algebra of 
the gauge group SU(N C ), while the Wilson loop operator Pexp (ig JdxA(x)) lies in its fundamental representation. 
The cusp anomalous dimension T cusp can be found from the multiplicative renormalization of the Wilson integral (7) 
[15, 37]: 

W(C x :iJ> 2 /\ 2 ,a s (i2 2 )) = Z cusp {C x -fl[i 2 ^ s {[i 2 )) ■ W bare (C x : fl 2 / X 2 , a s ) , (8) 

where ji 2 is the UV cutoff, fi 2 is the normalization point, and A 2 is the IR cutoff. The presence of the IR divergence in 
(8) is a common feature of on-shell amplitudes in massless QCD. Since Wb are knows nothing about the normalization 
point (the latter is fixed by choosing a concrete Z cusp ), one can write: 

\ r cU s P (C x ;a s (n 2 )) = -^ 2 ^\nW{C x :^ 2 /X 2 , as {p 2 )) = - f i 2 ^lnZ cusp (C x ;fi 2 /fi 2 ,a s (^)) . (9) 
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It can be shown that the cusp anomalous dimension (6) is linear in \ to all orders of perturbation theory in the 
large-g 2 regime [15]: 



r cusp (C x ;a s ) = lng 2 T cusp (a s ) + 0(ln° q 2 ) . 
Then, from the Eqs. (6, 9, 10) one finds after the simple calculations [14]: 

d\nF H {q 2 ) f^d^ 



d In q 2 



2£ 



1 cusp 



(10) 



(11) 



d\nF s (q 2 ) 



d In q 2 

where the "integration constant" of the hard part reads 



J x2 ^-pMOR 



d\nW np (q 2 ) 
d In q 2 



n<*.tf)) = 



and W np arises as the initial value of the soft part: 



2 .. _ d\nF H (q 2 ) 



9 In q 2 



d\nW np {q 2 ) _ dliiF s (q 2 



9 In q 2 



d In q 2 



(12) 



(13) 



(14) 



and is the only quantity where, according to our suggestion, the nonperturbative effects could take place [39, 40]. 
Then we get the Revolution equation of the total form factor at large q 2 : 



In 



F q (q 2 ) 



f 2 ^ 

hi 2 £ 



ln|- r coap (a a (0)-2r(a a (0) 



ln ^L^ Was(0)+ln ^iy 



• (15) 



In the next Section, we calculate explicitly the perturbative quantities entering Eq. (15) in one loop approximation. 



III. ANALYSIS OF THE PERTURBATIVE CONTRIBUTIONS TO THE WILSON INTEGRAL 

The analysis of the hard contributions [13, 14] at large q 2 yields: 



d\nF H (q 2 /ii 2 ,a s ) _ a 



d In q 2 



where Cf = (N 2 — l)/2N c . For the hard "integration constant" (13) one has: 



r(«.fe»)) = §^ 

4 7T 



(16) 



(17) 



The expression (16) is IR-safe, while the low-energy information is accumulated in the soft part of the quark form 
factor Fs- The Wilson integral (7) can be presented as a series: 



c n _2 JC X JC X JC X 

■ 6(x n ,x n -\...,x 1 ) Tr \A lin {x n )A lin _ 1 {x n - 1 )...A lil (x 1 )] |0> 



(18) 



where the function 9{x) orders the color matrices along the integration contour. In the present work, we restrict 
ourselves with the study of the leading order (one-loop — for the perturbative gauge field and weak-field limit for the 
instanton) terms of the expansion (18) which are given by the expression: 



w£l e (C x ) = -^ f dxj dy v D^(x-y), 



(19) 
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where the gauge field propagator D^ v (z) in n-dimensional space-time (n = 4 — 2s) can be presented in the form: 

D^(z)=g^d 2 A 1 (s,z 2 ,p 2 /X 2 )-d^d v A 2 (s,z 2 ,jl 2 /X 2 ) . (20) 

Here is the parameter of the dimensional regularization. The exponentiation theorem for non-Abclian path-ordered 
Wilson integrals [17] allows us to express (to one-loop accuracy) the Wilson integral (7) as the exponentiated one-loop 
term of the series (18): 

T4W(C x ;e,M 2 A 2 ) = cxp [\V^ re (C x ; e, fi 2 /X 2 ) + 0(a 2 s )] . (21) 

In general, the expression (19) contains ultraviolet (UV) and IR divergences, that can be multiplicatively renormalized 
in a consistent way [37]. In the present work, we use the dimensional regularization for the UV singularities, and 
define the "gluon mass" A 2 as the IR regulator. The dimensionally regulated formula for the leading order (LO) term 
(19) can be written as [39]: 

W^l e (C x ,s,fi 2 /X 2 ,a s ) = 8wa s C F h( X )(l- s)A 1 (s,0,fl 2 /X 2 ) , (22) 
where h(x) is the universal cusp factor: 

h( X ) = xcothx - 1 (23) 

and for perturbative gauge field 

4 *.ww=-^(4)'a. (24) 

The independence of the expression (22) from the function A 2 is a direct consequence of the gauge invariance. Then, 
in the one-loop approximation, 

W bare (C x ;e,p 2 /X 2 ,a s ) = l-^C F h( X ) Q - 1e + ln47r + In ^ , (25) 

and the cusp dependent renormalization constant, within the MS-scheme which fixes the UV normalization point, 
reads: 

Z cusp (C x ;e,fl 2 / f i 2 ,a s (^)) = l + ^lc F h(x) Q - lE + In 4^ + 0{a 2 s ) . (26) 

Using the Eq. (22), one finds the known one- loop result for the perturbative field, which contains the dependence 
on the UV normalization point /z 2 and IR cutoff A 2 (e.g., [15, 19]): 

W#> (C x ; ^ 2 /A 2 , a s (p 2 )) = -^lc F h( X ) In ^ + 0(a 2 s ) . (27) 

Therefore, in the leading order the kinematic dependence of the expression (19) is factorized into the function h(x), 
which at large-g 2 is approximated by: 

Mx)cxln^. (28) 

From the one- loop result (27), the cusp anomalous dimension which satisfies the RG equation (9) in one- loop order 
is given by: 

r£ sp (a s (p 2 )) = ^C F . (29) 

Substituting the anomalous dimension (29) in the one-loop approximation for the strong coupling into the Eq. (15), 
one finds 



exp 



2C F 



ln « 2 ,\ 3. lng 2 „ / InggN- 



+ W np (q 



F^(q 2 ). (30) 



In A 2 7 2^\nq 2 ' ^ w \~ ^lnA 2 , 

Note, that the exponent in Eq. (30) has an unphysical singularity at A 2 = 1 (in dimensional notations, A 2 = Aq CD ), 
i. e., where the coupling constant a s (A 2 ) has the Landau pole. This feature can be treated in terms of the IR 
renormalon ambiguities (see Section V), and is considered often as a signal of nonperturbative physics. In the present 
paper, we will consistently separate the sources of nonperturbative effects which can be attributed to uncertainties 
of the perturbative series resummation, from the "true" nonperturbative phenomena. An important example of the 
latter is provided by the instanton induced effects within the ILM of QCD vacuum, which is considered in the next 
Section. 
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IV. LARGE-g 2 BEHAVIOUR OF THE INSTANTON INDUCED CONTRIBUTION 



The instanton induced effects in the high energy QCD processes had been studied actively since the seventies 
[42, 43]). Recently the investigation of these effects renewed from the fresh points of view [34, 38-41, 44-47]. The 
Wilson integral formalism is considered as a useful and convenient tool in the instanton calculations, mainly due to the 
significant simplification of the path integral evaluation for an explicitly known gauge field [44]. Another important 
feature of this approach is the possibility of a correct analytical continuation of the results obtained in the Euclidean 
picture (where the instantons are only determined) to the realistic Minkowski space-time, where the scattering process 
(2) actually takes place [48]. Namely, in the instanton calculations, one maps the scattering angle, \i to the Euclidean 
space by the analytical continuation 

X^7, (31) 

and performs the inverse transition to the Minkowski space-time in the final expressions in order to restore the (in- 
dependence. Let us consider the instanton induced contribution to the function W np (q 2 ) from Eq.(14). The instanton 
field is given by 

A^x; p) = A«(x; p)^- = -R ab a a V ±b ^(x - z a ) v ^{x - z ; p), (32) 
z 9 

where R afc is the colour orientation matrix (a = 1, (TV 2 — l),b= 1,2,3) which provides an embedding of SU(2) 
instanton field into SU(3) colour group, a a, s are the Pauli matrices, and (±) corresponds to the instanton, or anti- 
instanton. The averaging of the Wilson operator over the nonpcrturbative vacuum is performed by the integration over 
the coordinate of the instanton center z , the color orientation and the instanton size p. The measure for the averaging 
over the instanton ensemble reads dl = dTl d^zo dn{p) 1 where dR refers to the averaging over color orientation, and 
dn(p) depends on the choice of the instanton size distribution. Taking into account (32), we write the Wilson integral 
(7) in the single instanton approximation in the form: 

wi(C 7 ) = ^-<0|TrexpKV a ) |0) , (33) 

where 



f) a (z , p) = R afc ?7 ± p,y / (x - z ) v (p(x - z ; p) . 

Jc~, 



(34) 



We omit the path ordering operator V in (33) because the instanton field (32) is a hedgehog in color space, and so it 
locks the color orientation by space coordinates. One obtains the all-order single instanton contribution in the form 
[39]: 

w J {C 1 )= Jd^zo Jdn(p) [cos ^(7, z , p) - cos 0(0, z , p)] , (35) 
where the squared phase </> 2 = <p a <p a may be written as 

4> 2 {i,z 0l p) = 

= V [{v lVj )zl - (v iZo )(v jZo )] / da V [((-iy +1 v % a-z ) 2 -p\ / da' V +1 v, a 1 

i,i=i,2 Ja Jo 



z ) 2 ;p] , (36) 



where — Pi^/m, and v\ 2 ~\, {v\V2) = COS7 in Euclidean geometry. Although sometimes such integrals (Eq. 
(36)) can be evaluated explicitly, the given contour requires numerical calculations, thus we restrict ourselves with the 
weak-field approximation which can be studied analytically. Then, in case of the instanton field, the LO contribution 
in Minkowski space reads 

wf\C x ) = 2h( X ) J dn{p) A{(0, p 2 A 2 ) , (37) 



where 



A((O,p 2 A 2 ) = -| I I e^^ 2 )[2^(fc 2 ;p)] 2 , (38) 
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and we use the same IR cutoff A 2 , while the UV divergences do not appear at all due to the finite instanton size. 
Here, ip(k 2 ;p) is the Fourier transform of the instanton profile function ip(z 2 ;p) and ip'(k 2 ;p) is it's derivative with 
respect to k 2 . In the singular gauge, when the profile function is 



V {z 2 ) = 



z 2 {z 2 + p 2 ) 



one gets: 



where 



A((0,p 2 A 2 ) 



^[ H p 2 \ 2 )^{p 2 X 2 ) + ^{p 2 \ 2 )\ , 



2 \2\ 



WA 



dz 



and 



$i(p 2 A : 



n=l J ° 



p 4 A 4 J z(l - z) 

[-p 2 \ 2 {xz + y(l-z))] 



2e 



2 \ 2 

z-p A 



lim $ (p A ) = 1 

A 2 ^0 



dxdydz 



n'.n 



_ c p 2 \ 2 [xz+y(l-z)] Um $1 ( p 2 A2) = q 

A 2 ->0 



are the IR-finite expressions. At high energy the instanton induced contribution is reduced to the form: 
d\nWj{q 2 ) _ 7T 2 



d In q 2 



y J dn(p) p 4 [ln(p 2 A 2 ) $ (p 2 A 2 ) + d^A 2 )] = -i? 7 (A 2 ) . 



Here we used the exponentiation of the single-instanton result in a dilute instanton ensemble [39] : 

Wi = exp (wi) , 



(39) 



(40) 



(41) 



(42) 



(43) 



(44) 



and took only the LO term of the weak-field expansion (19): = wj + {higher order terms). 

In order to estimate the magnitude of the instanton induced effect we consider the standard distribution function 
[50] supplied with the exponential suppressing factor, what has been suggested in [51] (and discussed in [52] in the 
framework of constrained instanton model) in order to describe the lattice data [53]: 



dn(p) = ^ C Nc 



2tt 



a s {n r )_ 



2N C 



exp 



2tt 



a s (p r ) 



(ppr) 13 exp (— 2nap 2 ) 



(45) 



where the constant C Nc = 4.6/tt 2 exp(-1.679iV c )/ [{N c - l)\(N c - 2)!] w 0.0015, a is the string tension, [3 = (3 + 
0(a s (p. r )), and p r is the normalization point [54]. Given the distribution (45) the main parameters of the instanton 
liquid model — the instanton density n and the mean instanton size p — will read: 



Jo 



dn(p) 



C N T((3/2 - 2) 



2tt 



2N C 



A 



QCD 



(2ira) 2 



(46) 



/ °°p dn(p) _ T(/3/2 - 3/2) 1 



f °°dn(p) T(P/2-2) 



(47) 



In Eq. (47) we choose, for convenience, the normalization scale p r of order of the instanton inverse mean size p^ 1 , 
taking into account that the distribution function (45) in the RG-invariant quantity up to 0(a 2 ) terms [54]. Note, 
that these quantities correspond to the mean size p and density n of instantons used in the model [28] , where the 
size distribution (45) is approximated by the delta-function: dn(p) = n S(p — po)dp . 
Thus, we find the leading instanton contribution (43) in the form: 



T s 9--4i 27T(7 

Ktt np In — — 



where 



K = 



r(/3 /2)[r(/3 /2-2)] ; 
2 [r(/3 /2 - 3/2)] 4 



0.74 



(48) 



(49) 
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FIG. 1: The asymptotic behaviour of the quark form factor is shown as the function of the dimensionless variable q 2 — Q 2 /A 2 , 
up to terms 0(fnlng 2 ). The long-dashed (red) presents the contribution of one loop perturbative terms; the solid (green) line 
represents the total form factor including the instanton induced part, Eq. (53). For comparison, the leading (~ In q 2 In In q 2 ) 
perturbative contribution is shown separately — the short-dashed (blue) line. 



and we used the one loop expression for the running coupling constant 

2tt „ HAL 



a s {p~ 



A) In pA 



2n t 



(50) 



The packing fraction TT 2 np 4 characterizes diluteness of the instanton liquid and within the conventional picture its 
value is estimated to be 0.12 , if one takes the model parameters as (see [29]): 



1/m 4 , p 1/3/m 



(0.44 GeVf 



(51) 



The leading logarithmic contribution to the quark form factor at asymptotically large q 2 is provided by the (pertur- 
bative) evolution governed by the cusp anomalous dimension (29). Thus, the instantons yield the subleading effects 
to the large-g 2 behaviour accompanied by a numerically small factor as compared to the perturbative term: 



Bi w 0.02 « 



2C F 



0.24 . 



(52) 



Therefore, from Eqs. (43) and (62), we find the expression for the quark form factor at large-g 2 with the one-loop 
perturbative contribution and the nonpcrturbative contributions (the function W np in Eq. (30)) which include both 
the instanton induced terms: 



F q (q 2 ) = exp 



2Cf . 2 , , 2 
In q In In q 

Po 



W [Bi 



2Cf\ 
Po J 



Coining 2 ) 



Fo(qh A 2 



(53) 



It is clear, that while the asymptotic "double-logarithmic" behaviour is controlled by the perturbative cusp anomalous 
dimension, the leading nonperturbative corrections result in a finite renormalization of the subleading perturbative 
term (Fig.l). Note, that the instanton correction has the opposite sign compared to the perturbative logarithmic 
term. 



V. AMBIGUITIES OF THE PERTURBATIVE RESULT: IR RENORMALONS AND ANALYTIZATION 

OF THE COUPLING CONSTANT 

As it was pointed out in the end of Section II, the perturbative evolution equation (30) possesses an unphysical 
singularity at the point A 2 = 1. Therefore, it is instructive to study the consequences of this feature. It is known that 
the presence of the Landau pole in the (one-loop) expression for the coupling constant leads to the IR renormalon 
poles in the Borel plane [21] which result in the renormalon-induced power corrections. The latter, being beyond the 
perturbative evaluations, are treated as a signal of nonperturbative effects which should compensate the ambiguities 
arising due to these poles. In the present situation, one can expect the corrections proportional to the powers of both 
scales: p? and A 2 . We assume here that the power /z 2 -terms are too strongly suppressed in large-g 2 regime and thus can 
be neglected in the given context, and focus on the power A 2 -corrections. To find them, let us consider the perturbative 
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function Ai(e, 0, jl 2 /\ 2 ) in the Eq. (22). The insertion of the fermion bubble 1-chain to the one-loop order expression 
(19) is equivalent to replacement of the frozen coupling constant g 2 by the running one g 2 — ► g 2 (k 2 ) — 47ra s (/c 2 ) [22] 
(for convenience, we work here in Euclidean space): 



A 1 (e,0,A2 2 A 2 ) = -4^ : 



2e 



-a s (k z ) 



ikz 



5{z 2 ) 



(2tt)™" sv " J k 2 (k 2 + X 2 ) ' 

Using the integral representation for the one-loop running coupling a s (k 2 ) = J^°da(l/k 2 )' Tb , b — /3o/47r, we find: 



(54) 



A 1 (e,0,f/X 2 ) = - 



A)(l - <0 



00 , HI - x - e)T(l + x + e) 

ax 



(x + e)T{l-e) 



(55) 



To define properly the integral in r. h. s. of Eq.(55), one needs to specify a prescription to go around the poles, which 
are at the points x n = n, n G N. Thus, the result of integration depends on this prescription giving an ambiguity 
proportional to (l/A 2 ) for each pole. Then, the IR renormalons produce the power corrections to the one-loop 
perturbative result, which we assume to exponentiate with the latter [22]. Extracting from (55) the UV singular part 
in vicinity of the origin x = 0, we divide the integration interval [0, oo] in two parts: [0,6] and [6, oo], where 6 < 1. 
This procedure allows us to evaluate separately the ultraviolet and the renormalon-induced pieces. For the ultraviolet 
piece, we apply the expansion of the integrand in Ai in powers of small x and replace the ratio of T- functions by 
exp(-7 E e): 



ArM,MW = - 



1 „ (ln47r- 7B + ln£r 

— t E (-) - 



/3o(l - e) 



kle n - k+1 



k,n=0 

which after subtraction of the poles in the M5-scheme becomes: 

1 ^ / y 2 X « ( _j„ 



Sl dxxn (h 



Ar(0,M7A 2 )- 



n=l 



».! 



dxx 



n-l 



In analogy with results of [49] , this expression may be rewritten in a closed form as 



Ar(0, M W 



1 



dx 



/3 (1 — e) Jo x 



-x In fi 



-a; In A 



Then, using the relation 



one finds 



dfi 2 



«"5 



d In q 2 

d\nW^\q 2 ) 
d In q 2 



2C F (l-e)Ar(0,M7A 2 ) , 



-irW p (a s ( M 2 ))(\-exp 



47T 



f3oce s (n 2 



(56) 



(57) 



(58) 



(59) 



(60) 



The second exponent in the last equation yields the power suppressed terms (l/g 2 ) in large-q 2 regime. In the leading 
logarithmic approximation (LLA) Eq. (59) is reduced to: 



dW^jq 2 ) 
din q 2 



2cy 



In 



In /i 2 

LTa 2 



(61) 



The last expression obviously satisfies the perturbative evolution equation (30). 

The remaining integral in Eq. (55) over the interval [6, oo] is evaluated at e = since there are no UV singularities. 
The resulting expression does not depend on the normalization point /i, and thus it is determined by the IR region 
including nonperturbativc effects. It contains the renormalon ambiguities due to different prescriptions in going 
around the poles x n in the Borcl plane which yields the power corrections to the quark form factor. 

After the substitution [i 2 = q 2 and integration, we find in LLA (for comparison, see Eq. (30)): 



jpren 



(<1 2 



cxp 



-^hV (lnln 9 2 -l)-hV$ re „(A 2 ) 



(62) 
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where the function $ ren (A 2 ) = X)fe=o ( t ) k{^/^ 2 ) k accumulates the effects of the IR renormalons. The coefficients 
cannot be calculated in perturbation theory and are treated often as "the minimal set" of nonpcrturbative parameters. 
It is worth noting that the logarithmic g 2 -dependence of the renormalon induced corrections in the large-q 2 regime is 
factorized, and thus the Eq. (62) corresponds to the structure of nonperturbative contributions found in the one-loop 
evolution equation (30), in a sense of its largc-g 2 behaviour. On the other hand, the IR structures of the renormalon 
corrections and the instanton induced ones (43) are different. In our point of view, it allows us to separate the true 
nonperturbative {e.g., instanton induced, but not only) effects from that ones related to ambiguities of the resumed 
perturbative series. 

The latter conclusion can be illustrated by considering the consequences of an analytization of the strong coupling 
constant [23] in the perturbative evolution equation. In this approach, the one- loop strong coupling a s (/i 2 ) is replaced 
by the expression which is analytical at p? = 1 (i.e., at Aqcd in dimensional variables): 



..AJV/..2N 4tt / 1 . 1 



The direct substitution of (63) into the evolution equation (15) yields (for brevity, we assume q 2 = A 2 ): 



2C F 



\nF q AN (q 2 ) = 



2 , In q 2 t 9 3 , In q 2 



= In q In - — — — In — — — In ■ 



+ In q 2 ( In * + In *±) \ (in 2 q 2 - In 2 A 2 ) 



In A 2 A 2 2 in A 2 
q 2 A 2 -l\ 1 

-Li 2 (l - q 2 ) + Li 2 (l - A 2 ) + | (in + In . (64) 

The functions Li 2 in the resulting expression accumulate the power corrections of q 2 as well as IR scale A 2 , but does 
not exhibit a singularity at A 2 = 1. Therefore, it gives no room for IR renormalons ambiguities, at least in the 
considered approximation. Nevertheless, the power corrections of a nonperturbative origin do contribute to the large- 
q 2 behaviour, and the investigation of the correspondence between latter and the instanton corrections calculated in 
the previous Section would be an interesting task. Note, that the consequences of the analytization of the strong 
coupling constant in the IR region have been studied earlier in the case of the Sudakov effects in the pion form factor 
and Drell-Yan cross section in the works [55]. 

Another possible way to avoid the Landau pole on the integration path have been developed within the dimensional 
regularization [3]. In this case, the running coupling reads 

"■^'MW^V (65) 

and for complex e , Re e < it has the Landau pole at the complex value of /x 2 , that is this singularity appears to 
be out of the integration contour. In the limit e — > 0, the form factor reads [3] (for comparison, see Eq. (30): 



FP R tf) = exp 



.XjL ^ + In e (in q 2 - + In q 2 (in In q 2 - l) - | In In q 2 ^j + 0(e, e In e) 



(66) 



This expression also leaves no room for any renormalon induced effects. In the same time, the instanton induced 
contribution still takes place since they enter into the "integration constant" W np which knows nothing about the 
analytical properties of the coupling. It should be emphasized that the absence of the IR renormalon induced power 
corrections obviously does not mean the absence of the power corrections at all. It merely implies that the relations 
between the ambiguities in the perturbative resummation procedures and the true nonperturbative physics are not so 
evident and should be studied in more detail. 



VI. CONCLUSION 



The structure of the nonperturbative corrections to the quark form factor at large momentum transfer was analyzed. 
In this work, the quarks were assumed to be on-mass shell. In order to model the nonperturbative effects, the quark 
scattering process was considered in the background of the instanton vacuum. The instanton induced contribution 



11 



to the electromagnetic quark form factor is calculated in the large momentum transfer regime. It was shown that 
the instanton induced corrections correspond to the leading term proportional to hiq 2 . The magnitude of these 
corrections is determined by the small instanton liquid packing fraction parameter, and they can be treated as finite 
renormalization of the subleading logarithmic perturbative part (53). 

We have to comment that the weak-field limit used in the instanton calculations may deviate from the exact 
result. Nevertheless, we expect that using of the instanton solution in the singular gauge, that concentrate the field 
at small distances, leads to the reasonable numerical estimate of the full effect. Thus, the resulting diminishing of 
the instanton contributions with respect to the perturbative result appears to be reasonable output. It should be 
emphasized that in the present paper, all the calculations have been performed analytically while the evaluation of 
the instanton contribution beyond the weak-field approximation requires the numerical analysis, what will be the 
subject of a separate work. Besides this, the use of the singular gauge for the instanton solution allows one to prove 
the exponentiation theorem for the Wilson loop in the instanton field [39] which permits to express the full instanton 
contribution as the exponent of the all-order single instanton result (44). 

It is also important to note that the results are quite sensitive to the way one makes the integration over instanton 
sizes finite. For example, if one used the sharp cutoff then the instanton would produce strongly suppressed power 
corrections like oc (l/g) /3 °. However, we think that the distribution function (45) should be considered as more 
realistic, since it reflects more properly the structure of the instanton ensemble modeling the QCD vacuum. Indeed, 
this shape of distribution was recently advocated in [51, 52] (see also [31] and references therein) and supported by 
the lattice calculations [53] (for comparison, see, however, [56, 57]). 

The instanton induced effects are more interesting for investigation and more important for phenomenology in the 
hadronic processes which possess two energy scales, such as the total center-of-mass energy s (hard characteristic 
scale), and the squared momentum transfer —t which is small compared to the latter: —t « s. One of the most 
interesting examples of such processes is the parton-parton scattering and the soft Pomeron problem [4, 44, 45]. 
Another important situations where the nonperturbative (including instanton induced) effects can emerge are the 
transverse momentum distribution of vector bosons in the Drell-Yan process (see, e.g., [22]), and the phenomenon of 
saturation in deep-inelastic scattering (DIS) at small-x [58]. The explicit evaluation of the instanton effects in some 
of these processes will be the subject of our forthcoming study. 
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